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ABSTRACT 

We numerically compute light and polarisation curves of 7-ray burst afterglows for 
various configurations of the jet luminosity structure and for different dynamical evolu- 
tions. We especially consider the standard homogeneous "top hat" jet and the "univer- 
sal structured jet" with power-law wings. We also investigate a possible more physical 
variation of the "top hat" model: the "Gaussian jet" . The polarisation curves for the 
last two jet types are shown here for the first time together with the computation 
of X-ray and radio polarised fluxes. We show that the lightcurves of the total flux 
from these configurations are very similar to each other, and therefore only very high 
quality data could allow us to pin down the underlying jet structure. We demonstrate 
instead that polarisation curves are a powerful means to solve the jet structure, since 
the predicted behaviour of polarisation and its position angle at times around the jet 
break are very different if not opposite. We conclude that the afterglow polarisation 
measurements provide clear footprints of any outflow energy distribution (unlike the 
lightcurves of the total flux) and the joint analysis of the total and polarised flux 
should reveal GRBs jet structure. 

Key words: gamma-ray: bursts — radiation mechanisms: non thermal — polarisa- 
tion 



1 INTRODUCTION 

After the discovery of a significant, albeit small, degree 
of linear polarisation of the afterglow optical flux of GRB 
9905f0 (Covino et al. f999; Wijers et al. 1999), there have 
been several other detections of linear polarisation in a num- 
ber of afterglows (see Covino et al. 2002 for review). Typ- 
ically, the polarisation is observed to be at the 1-3 per- 
cent level (but see Bersier et al. 2003 for GRB 020405), 
with constant or smoothly variable level and position angle 
when associated with a relatively smooth lightcurve (e.g. 
GRB 020813, Gorosabel et al. 2004). When deviations from 
a smooth power-law decay in the lightcurve are instead 
present, polarisation curves show a certain degree of com- 
plexity (e.g. GRB 021004 Lazzati et al. 2003; Rol et al. 2003; 
GRB 030329 Greiner et al. 2003). 

In order to observe polarisation, some kind of asymme- 
try is needed: this can be provided by patches of coherent 
magnetic field, as suggested by Gruzinov & Waxman (1999) 
and Gruzinov (1999). In addition, small regions in which 
the magnetic field has some degree of order could be ampli- 
fied by scintillation (Medvedev & Loeb 1999), or by gravita- 
tional mini-lensing (Loeb & Perna 1998; Ioka & Nakamura 



2001). In these cases the required degree of asymmetry is 
in the structure of the magnetic field, and not in the over- 
all geometry of the fireball, which could even be spherically 
symmetric. 

Recently Granot & Konigl (2003) proposed that the 
required asymmetry might be provided by an ordered mag- 
netic field embedded in the circum-burst material, possibly 
amplified when the shock propagates in it, to reach values 
close to equipartition with the energy density of the shocked 
material in the "pulsar wind bubble" scenario (Konigl & 
Granot 2002). 

A different class of models postulates that the fireball 
is jetted (i.e. the ejecta are collimated into a cone with half 
opening angle 6 jet). In this case the observer likely sees the 
fireball off-axis, since the probability to be exactly on-axis 
is vanishingly small (Ghisellini & Lazzati 1999; Sari 1999, 
hereafter GL99 and S99, respectively). When the fireball 
bulk Lorentz factor r is ~ l/(#jet — 8 ) (where o is the 
viewing angle) the emitting surface starts to be asymmetri- 
cal with respect to the line of sight. Moreover it is assumed 
that a magnetic field is compressed in the plane normal to 
the motion, analogous to what has been proposed for the 
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Figure 1. Cartoon of the three jet configurations discussed in 
this paper. The figure shows the energy per unit solid angle of the 
jets logaritmically scaled. The scale is different for each jet: it has 
been chosen in order to visually emphasize the characteristics of 
each configuration. 



jets of radio-loud active galactic nuclei (Laing 1980). Face- 
on observers would see a completely tangled magnetic field, 
but edge-on observers (in the comoving frame) would see 
an aligned field, and therefore would detect synchrotron po- 
larised radiation. Since the fireball is moving with a high 
bulk Lorentz factor, the edge-on comoving observer corre- 
sponds to an observer in the lab located at an angle ~ 1/r. 
According to this idea, there is a tight link between the be- 
haviour of the total and the polarised flux as a function of 
time. The light curve of the polarised flux has two maxima 
(corresponding to the presence of the two edges of the jet) 
with a position angle switched by 90°. The maxima occur 
just before and after the achromatic "jet" break in the light 
curve of the total flux. A main assumption of this model is 
that, at any given angle from the apex of the jet, the lu- 
minosity emitted per unit solid angle along the jet axis and 
along the jet borders is the same. Let us call jets with this 
energy structure "homogeneous" jets (HJs) (see Fig.l). 

It is possible, instead, that the radiated power (per unit 
solid angle) along the jet axis is larger than what is emitted 
along the "wings" . If the wing energy distribution is a power- 
law, we refer to these configurations as "structured" jets 
(SJs) (see Fig.l). As Rossi, Lazzati & Rees (2002) (thereafter 
RLR02) and Zhang & Meszaros (2002) have demonstrated, 
if the luminosity per unit solid angle is L{9) cx 6~~ a with 
a close to 2, then observers with a viewing angle o would 
see an achromatic jet break when F ~ l/9 . It is therefore 
possible that all GRB jets are intrinsically alike, having the 
same total intrinsic power and the same jet aperture angle: 
they appear different only because they are viewed along 
different orientations. If the jet were uniform, instead, they 
should have a large variety of aperture angles, to account 
for the different observed jet-break times (Frail et al. 2001). 



As demonstrated analytically in RLR02 and more re- 
cently numerically by Salmonson (2003) (thereafter S03), 
it is difficult, on the basis of the observed light curve, to 
discriminate between homogeneous and structured jets (see 
also Granot & Kumar 2003). However, as described in this 
paper, the two models are markedly different in the polar- 
isation properties of the produced afterglow flux. In both 
models the polarisation is produced because different parts 
of the emitting jet surfaces do not contribute equally to the 
observed flux. In the homogeneous jet model this starts to 
occur when 1/r becomes of the order of 6 jet — 9 (i.e. when 
the emitting surface available to the observer "touches" the 
near border of the jet). In the structured jet model, instead, 
the required asymmetry is built-in in the assumption that 
the emission is a function of 6, so that the relevant emitting 
surface is never completely symmetric for off-axis observers. 

Finally we consider a jet with a Gaussian luminosity dis- 
tribution (Zhang & Meszaros 2002). This can be regarded 
as a more realistic version of the sharp edged standard jet: 
the emission drops exponentially outside the typical angu- 
lar size (6 C ), within which it is roughly constant (see Fig.l). 
Let us call it the "Gaussian jet" (GJ). It has been argued 
that this configuration can accomodate a unified picture of 
GRBs and X-ray flashes. The underlying assumptions is the 
presence of an emission mechanism for which the peak en- 
ergy E p in the prompt emission is a decreasing function of 
the angular distance from the jet axis. In this way X-ray 
flashes would be the result of observing a GRBs jet at large 
angles. According to Zhang et al. (2003) the GJ would re- 
produce the observed correlation E p cx E]^ (Amati et al., 
2002), while under the same assumption the universal SJ 
and the "top hat" jet would face severe problems. We no- 
tice here, however, that the above correlation is still based 
on a very small database in the X-ray flashes regime (only 
two X-ray flashes are included) and should be confirmed by 
future observations. As regards afterglow properties, a GJ 
seen within the core produces lightcurves that are similar 
(but with smoother breaks) to the HJ's ones (Granot & Ku- 
mar 2003). The luminosity variation with angle gives, as in 
the case of a SJ, a net polarisation without the need of edges 
and we show here that its temporal behaviour is indeed dif- 
ferent from both the SJ's and the HJ's one. 

The detailed analysis of the polarisation characteristics 
and their connection with lightcurves in the homogeneous 
jet, universal structured jet and in the Gaussian jet models 
is the main goal of this paper. In Section 2 we described the 
numerical code we have implemented to study these mod- 
els. Some analytical and semi-analytical results have been 
derived by GL99 and S99, for a non-spreading jet and for a 
sideway expanding jet respectively. The simplified prescrip- 
tion assumed by S99 for the lateral expansion led to predic- 
tion of a third peak in the light curve of the polarised flux 
for an observer close to the border of the jet. In addition, 
GL99 did not consider, for simplicity, the effects of the dif- 
ferent travel times of photons produced in different regions 
of the fireball, while S99 considered this effect by represent- 
ing the viewable region as a thin ring centred around the 
line of sight: the ring has an angular size of and a 

constant width with respect to the ring radius. Our numer- 
ical approach allows us to include the effects of the different 
photon travel time and to analyze and compare different 
prescriptions for the side expansion of the fireball. 
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This paper is organised as follows. In §2 we present the 
model for the jet and for the magnetic field. In §3 we show 
the results for a homogeneous jet, in §4 those for a structured 
jet and in §5 those for a Gaussian jet. The comparison and 
discussion can be found in §6. Finally in §7 we derive and 
discuss our conclusions, adding possible complications to the 
models. 

Throughout this paper the adopted cosmological parameters 
are Ho = 65, fl\ = 0.7 and Q m = 0.3. 



2 THE CODE 

2.1 The Jet structure and dynamics 

In this paper we show results obtained with two different 
codes. The first one is fully described in S03 while the second 
one is discussed in this section. The main difference between 
the two codes is in the treatment of the sideway expansion 
and dynamics in the non-relativistic phase. In the follow- 
ing we only remind the reader of the different assumptions 
adopted by S03 and refer the reader to the paper for more 
details. 

We assume that the energy released from the engine 
is in the form of two opposite jets. They are described by 
the following distributions of initial Lorentz factor To and 
energy per unit solid angle e with respect to the jet axis 
{9 = 0): 
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where 9j et is the jet opening angle, 6 C is the core angular 
size, e c = e(0) and r c = T(0). In the following, in order to 
make the comparison with the homogeneous jet easier, we 
will use preferentially the local isotropic equivalent energy, 
defined as Ei 3O (0) = 4ne(6). In Eqs.^and^a,,, or controls 
the shape of the energy and To distributions in the wings, 
while f3 e ,Pr controls the smoothness of the joint between 
the jet core and its wings. 

If a e ,Qr = 0, equations Q and |21 describe the standard 
top hat model with sharp edges (homogeneous jet). If the 
observer line of sight is located within the jet, the observer 
detects the GRB prompt phase and its GRB afterglow (GA) ; 
if the viewing angle o is larger then 9j e t, he observes what 
it is called an orphan afterglow, an afterglow not preceded 
by the prompt 7-ray emission. 

When a e > 0, the code describes a structured jet; in 
this case 9j e t is assumed to be always much larger then the 
observer angle. In fact we consider here a boundless jet (the 
end of the wings are so dim that are undetectable), in con- 
trast to the sharp edged homogeneous jet. If a e = 2 and 
(3 e — » 00, the structured jet is that described in RLR02, 
while S03 adopts a e = 2 and f3 t — 1. 

For the Gaussian jet we use instead 
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For simplicity we assume axial symmetry. Our initial 
Lorentz factor distribution satisfies 1/r < 6, therefore re- 
gions on the shock front with different Fq and energy are 



causally disconnected and they evolve independently un- 
til 1/r = 6. In the numerical simulations we assume for 
all models a costant initial Lorentz factor across the jet, 
with To = 10 4 . With this choice the lightcurves shown 
in this paper (for observed times ^ 15 min) are insen- 
sitive to the initial V distribution. As a matter of fact, 
for any Tq{6) ^ 50, the fireball deceleration starts earlier 
(t d ~ 245s(£53/no/ro,2) 1,/3 ) than the smallest time of the 
figures and afterwards the evolution follows the BM self sim- 
ilar solution and consequently the shown afterglow proper- 
ties are independent from the initial Lorentz factor. If rela- 
tivistic kinematic effects freeze out the lateral expansion or 
the pressure gradient prevents mixing, the different parts of 
the flow are virtually independent along their entire evolu- 
tion. We allow each point of spherical coordinates (r, 8, (f>) to 
evolve adiabatically and independently, as if it were part of 
a uniform jet with e = e(6), To = Fo(0) and semi-aperture 
angle 6. Therefore, if the mixing of matter is unimportant 
this treatment is correct at any time, otherwise it gives an 
approximate solution for 1/f > A9. Actually numerical hy- 
drodynamical simulations seem to suggest that e{9) does not 
vary appreciably with time until the non-relativistic phase 
sets in (Kumar & Granot 2003) thus supporting our numer- 
ical approach. 

The full set of equations that determine the dynamics 
of each patch of the jet is: 



r = 



y / i + 4r / + 4/2-i 
2/ 



(4) 



(e.g. Panaitescu & Kumar 2000, thereafter PK00) where the 
parameter / (the ratio of the swept-up mass to the initial 
fireball rest mass) is given by: 



/ = 



M (e) j„ 



r 0.(r) p(r) dr, 



(5) 



where Mo is the rest mass of the two (symmetric) jets, 
fi(r) = 4-7T (1 — cos#(r)) is their solid angle and p(r) is the 
ambient medium matter density. The evolution of the solid 
angle is described by: 



dO 

dr 



Cs{6,r) 
c(3Yr 



(6) 



For the comoving lateral velocity c s we tested three different 
recipes. First, we analyze a non-sideways expanding (NSE) 
jet with 



0. 



(7) 



then a sideways expanding (SE) one, either with a constant 
comoving sound speed (Rhoads 1999) 



(8) 



or with a more accurate treatment, which takes into account 
the behaviour of the sound speed as a function of the shock 
Lorentz factor: 



7(7-i)(r-i) 

1+ 7 (T-1) 



(9) 



(Huang et al. 2000), where 7 = 

In the code developed by S03 (see in particular S03 §2), 
it is assumed that both momentum and energy are conserved 
(Rhoads 1999): 
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To + / 
V/1 + 2T0/ + / 2 



(10) 



This affects mainly the temporal slope of the lightcurve in 
the non-relativistic regime, since the jet does not follow the 
Sedov- Taylor solution. For the sideways expansion prescrip- 
tion S03 assumes that the lateral kinetic energy of the shock, 
in its radially comoving frame, is a constant proportion of 
the radial kinetic energy: 



Rk = W± - l)/(7 - 1), 



(11) 



where Rk (using S03's notation) describes the shock effi- 
ciency to convert radial kinetic energy (7 — 1) in lateral 
kinetic energy (7^ — 1) and 



r = 77I 



(12) 



The jet dynamics in the trans- and non-relativistic phase 
is uncertain. If the two opposite jets do not merge (9j e t is 
always < 90°) when the dynamics becomes non-relativistic, 
the radial momentum should be conserved through all the 
evolution of the jet. This possibility depends strongly on the 
initial opening angle and on the assumed type of lateral ex- 
pansion. For example, using the lateral velocity assumed by 
S03, initially narrow jets will not merge, because the side- 
ways expansion, deriving its energy from the forward expan- 
sion, is soon exhausted. On the other hand, if, at late times, 
lateral spreading causes the jets to homogenize and become 
effectively spherically symmetric, then Eq. gjis correct. For 
these reasons and the presence in literature of both treat- 
ments we compare and show in this paper results from both 
Eq. BJand Ed-ITU! 

In Fig. [5] we show the opening angle of the jet vs. 
the shock radius for all sideways velocity prescriptions. The 
jet spreads more efficiently when the simpler prescription 
c s — c/\/3 is adopted while with Eq. llll the expansion stalls 



out at fairly small angles (~ 0j (0) + arctan - 



since the 



jet is radially decelerating. Eg. 1111 gives another interesting 
behaviour: the jet begins to laterally expand earlier than 
9j e t(r) ~ 1/r; this is because Rk = 0.01,0.1 correspond to 
very large and supersonic initial lateral expansions. A direct 
comparison between the two sonic expansions (dotted and 
dashed lines) shows that a variable sound velocity (Eq. 
gives a final 6j e t smaller than Eq. |H| Note however that 
this seems not to affect appreciably the resulting lightcurve 
(Fig. I). 



2.2 Spectrum and Luminosity 

In order to isolate the effects of the jet structure on the emis- 
sion and polarisation curves, we assume throughout this pa- 
per a constant density environment p(r) = m v n (m p being 
the proton mass and n the number density). Polarisation 
curves derived here should be therefore compared only to 
data from afterglow with smooth lightcurves (Lazzati et al. 
2003). Polarisation curves in a windy environment are sim- 
ilar to ISM ones, but characterised by a slower evolution 
(Lazzati et al. 2004). 

The comoving intensity I' is assumed to be given by the 
synchrotron process only (Granot & Sari 2002); we do not 
include inverse Compton emission, because (for the adopted 
fiducial parameters) it can modify the observed afterglow 
lightcurve and polarisation curves only in the X-ray band 




r(cm) 

Figure 2. The opening angle of an homogeneous jet as a function 
of the distance from the explosion site for four different lateral 
velocities given by Eq. Eq. |S] Eq. |5]and Eq. [lj] The jet has 

an initial 6 Je t of 10°, To = 10 4 , Ei so = 10 53 and the external 
density is n = 10 cm -3 . See text for discussion. 



after ~ 10 days (depending on the density, see e.g. Sari & 
Esin 2001). Therefore our results are accurate up to the X- 
ray band. We will comment on how the polarisation curves 
would be modified, should IC emission be included (§3.3.4). 
In the shocked matter the injected electron Lorentz fac- 



tor is 



7m 



(r-i), 



(13) 



where m e is the electron mass and e e is the fraction of kinetic 
energy given to electrons; the Lorentz factor of the electrons 
that cool radiatively on a timescale comparable to dynamical 
timescale (i.e. time since the explosion) is 



7c 



15 tt- 



a T B>r < (14) 
where <tt is the Thomson cross section and the comoving 
magnetic field B is given by: 



B = ^32 7re B m p c 2 ny/r 2 - 1. 



(15) 



where e_g is the fraction of internal energy that goes to the 
magnetic field. Note that Eq. 1141 and Eq. 1151 are different 
from those given by PK00 that are not accurate in the trans 
and sub-relativistic regimes. The corresponding synchrotron 
frequencies are: 



0.25 



e^B 



(16) 



where the constant (0.25) is calculated for an electron en- 
ergy distribution index p = 2.5 (PK00), e is the electron 
charge, i — m for the peak frequency and i = c for the cool- 
ing frequency. The synchrotron self absorption frequency is 
given by: 

v a = Vi t? /o , (17) 

where i = m in the slow cooling regime and i — c in the fast 
cooling regime; Tj = ^ -^5- is the optical thickness at Vi. 

The comoving peak intensity I' p is 
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4 = 



P'nr 



(18) 



where P' = | ot c (7^ — l) is the total power emitted 
by relativistic electrons with Lorentz factor 7 m and isotropic 
distribution of pitch angles. The local observed luminosity 
is then computed through: 



dL(t, 8, cf>) = I 1 S 3 r 2 sin 



where 5 



dOd<j> (19) 
is the relativistic Doppler factor and 



r (1-/3 cos e) 

9 is the angular distance from the line of sight. The observed 
arrival time t of photons can be computed as follows: 



t — tlai — 



(20) 



where the time in the laboratory frame t; a j, at which the 
photons were emitted is 



tlab — 



dr 



(21) 



where /3 a h c is the shock front speed. The Lorentz factor T s h 
of this front is related to the Lorentz factor F of the shocked 
matter behind it (Eq.gl by F sh = 1 + a/I (T - 1) (e.g. Sari 
1997). 



2.3 The emitting volume 

For a given time, Eg. I2UI describes the locus of points from 
which photons arrive simultaneously at the detector (equal 
arrival time surfaces EATS). The EATS for a structured 
jet are shown in Fig. |3 where the observer is located to 
the right. Unlike the homogeneous case (e.g. Panaitescu & 
Meszaros 1998; Granot, Piran & Sari 1999), the EATS shape 
in the SJ depends on the viewing angle. This is because, for 
relativistic effects, each line of sight mimics a homogeneous 
jet with different parameters: 9j e t — 28 , Ei ao ~ Ei 30 (8 ) 
and To ~r (6» o ) (RLR02). 

The emission that the observer detects at time t does 
not come only from a thin layer near the EATS; in fact it 
comes from a sizable fraction of the volume behind it whose 
width is ~ — • Therefore, for each t, 9 and <f>, we integrate 
equation ll9l also over this emitting volume, using the Bland- 
ford & Mckee (1976, hereafter BM) solution for the Lorentz 
factor of the shocked gas behind the shock front (Granot, 
Piran & Sari 1999). This calculation for the observed flux is 
strictly valid for v > v a , since we do not take into account 
self absorbing effects. 



2.4 Magnetic field configuration and linear 
polarisation 

The magnetic field configuration we adopt is obtained by 
compressing, in one direction, a volume containing a ran- 
dom magnetic field: it has some degree of alignment seen 
edge-on while it is still completely tangled on small scales 
in the uncompressed plane. This could be the geometry of a 
magnetic field produced by the blastwave that, sweeping up 
the external medium, could confine the field in the sky plane. 
It may also be the natural configuration of shock generated 
magnetic fields (Medvedev & Loeb 1999). 

Since the fireball is relativistic, the circle centered 
around the line of sight, with angular aperture 9 ~ 1/r, 
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Figure 3. The Equal Arrival Time Surfaces for a a SJ seen at 
0° (solid line), 6° (dashed line) and 12° (dot-dashed line) from 
the jet axis. The straight dotted lines mark the position of the 
jet axis and the horizontal dashed line is the line of sight. For 
9 a = of course the dotted line overlaps the dashed line. The 
jet parameters are: external density n = 1 cm -3 , core dimension 



a jet 



0.1, E c = 10 54 erg, T c 



10 4 



2 and tip 



EATSs clearly depend on the observer position: each line of sight 
mimics an homogeneous jet with different energy, Lorentz factor 
and e jet . 



is the region contributing the most to the observed emission 
and the photons that reach the observer from the borders of 
that circle are emitted at 9' = 90° in the comoving frame. 
Therefore the light coming from an angle 9 ~ 1/F has the 
maximum degree of polarisation Po while the light travelling 
along the line of sight is unpolarised (see Fig. 1 in GL99). 
We assume Po = 60%. In the comoving frame this value 
decreases with the angular distance (9') to the line of sight 
as 



P{6) = Po 



1 + cos 2 



(22) 



(Laing 1980). The Lorentz transformations of angles give us 
the relation between 8 and 9': 



■(3 



1 - f3 cos e 
sin^ =rsin0 (l+0 cosfl') 



(23) 
(24) 



In terms of the Stokes parameters, Q, U, V, and local 
luminosity dL (equivalent to the Stokes parameter I) P(9) 
can be described as a vector with components: 



dQ(t,0 
dU{t,8 
and 
V = 



= P(9)dL(t,9,<f>) cos(2(j>), 
= P(9)dL(t,9,(j>) sin(20). 



(25) 
(26) 

(27) 



Integrating Eq.|2H]and Eg. 1251 over the EATS we obtain the 
intensity of the observed total polarisation vector at a time 
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Figure 4. Comparing the R band (7 X 10 14 iiz) lightcurves re- 
sulting from different sideways expansions shown in Fig.|2] where 
the jet parameters are given. The shock parameters are: p = 2.5, 
e e = 0.01 and eg = 0.005. An on-axis observer (0 O = 0) has been 
assumed here. SE jets with constant (dashed line) and variable 
(dotted-line) sonic lateral velocity give very similar lightcurves. 
SE jets with supersonic lateral velocities ( "dash-dot" and "dash- 
3 dots" lines) have been divided by a factor of 10 for clarity. 
They have more extreme spreading jet features (smoothness of 
the break, late time break) and due to a different dynamics they 
do not follow the Sedov- Taylor solution in the non-rclativistic 
regime. 

which is actually the fraction of flux linearly polarised. The 
direction of the total polarisation vector is then 

4>{t) = 0.5 arctan^(. (29) 



3 HOMOGENEOUS JET: RESULTS 

In Fig. Hand Fig. we compare the light and polarisation 
curves resulting from different recipes for the lateral expan- 
sion (§2.1). The expansions given by Eq. |S]and Eq. ^pro- 
duce very similar temporal behaviors; therefore in the fol- 
lowing we show results which are strictly valid for a variable 
sound speed but the same discussion holds for c s = c/v3. 
The lightcurves shown in Fig. 0] are all consistent with the 
same pre and post-break temporal slopes for a relativistic 
jet. However the curves corresponding to a jet that obeys 
the conservation of momentum (Eg. 1101 the dash-dotted and 
dash-3 dotted lines) does not follow a Sedov- Taylor model 

O/r (3p-4) 

(/? oc r ' and flux F(t) oc t 2 ( for v > u c ) in the non- 
relativistic regime but j3 oc r~ 3 and F(t) oc i -43 / 16 . Then in 
this latter case the flux falls off more rapidly while in the for- 
mer case the lightcurve tends to flatten. We anticipate that 
in both cases the non-relativistic slopes do not depend on 
the viewing angle but only on the spectrum, contrary to the 
temporal slopes in the relativistic regime (see also Fig.|S]|7|l. 
Fig. 2] and Fig. also show that a SE jet has 

• a smoother break in the lightcurve, 

• a later time break, 
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Figure 5. The same as Fig. [3] for the polarisation curves. SE 
jets with constant (dashed line) and variable (dotted-line) lateral 
velocity give almost indistinguishable polarisation curves (the sec- 
ond highest curve). Most of the time SE jets have a lower degree 
of polarisation than the NSE jet (solid line). Note that the larger 
the lateral velocity, the lower the polarisation. 

• smaller polarisation peaks, 

compared to a NSE jet and these characteristics are even 
more evident for jets with a supersonic sideways expansion. 

All these features (among others) are discussed in more 
details in the following, comparing a NSE jet with a jet 
undergoing a sideways expansion given by Eq. [§] For a jet 
evolution following Eq. llll we refer the reader to S03, where 
a more complete discussion is given. 

3.1 Lightcurve: the shape of the break 

The R band lightcurve for a NSE and for a SE homogeneous 
jet are shown in Fig. |S| and Fig. |7| respectively. They show 
interesting features, some of which have never been discussed 
before. In the lower panels the temporal index a (defined 
as F(t) oc t a ) is plotted versus time. We call ot\,ct2 and 
Q3 the pre-break, the post-break and the non-relativistic 
slope respectively. The horizontal dot-dashed lines show the 
expected slopes from on-axis standard calculations. 

3.1.1 a 2 

The "standard" post-break slopes (central dot-dashed lines 
in Fig.|Uand Fig.|7] lower panels) are calculated considering 
the loss of the emitting area from a spherical blast wave 
(1 — cosT -1 ) to a conical one (1 — cos jet)- For the SE 
jet we also consider the more rapid deceleration due to the 
increase of the shock front surface (Rhoads 1999; Sari, Piran 
& Halpern 1999). In both cases the surface brightness is 
supposed to be the same across the surface of the jet at all 
times. For v > v c the expected breaks are Aa nse = | for a 
NSE jet and Aa se = = 1.125 for a SE one. 

In fact we find that the flux after the break falls off 
more rapidly than expected, and this effect is more evident 
as the electron energy distribution N(-y) oc r y~ p becomes 
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Figure 6. Lightcurves (upper panel) and power-law indices as 
a function of time (lower panel) for a NSE homogeneous jet 
(a E = ap=0 and /3 e = /3p = 1), seen from different viewing an- 
gles. The parameters are To = 10 4 , Ei so = 10 53 erg, 9 jet = 10°, 
n = lOcm" 3 , e e = 0.01, e B = 0.005, p = 2.5, v = 7 X 10 14 Hz 
and 2 = 1. For these parameters the R band is beyond the cooling 
frequency. The curves corresponding to O = — + 9 = 0.89j e t 
are shown but not labelled one by one for clarity. The black 
solid line corresponds to 9 a = 9 jet, that shows a sort of in- 
termediate behavior between GRB afterglows (9 < 9 jet) and 
orphan afterglows (9 > 9j e t). The temporal decay for GRB 
afterglows is shown in the lower panel. The horizontal dash- 
dotted lines are the standard analytically predicted slopes for 
the pre- break power-law ct\ = — ^ 3p 4 2 ~* = —1.375, the post-break 
power-law c*2 = oil — | = 2.125 and the non-relativistic phase 
03 = — ( 3p 2 4 ) = 1.75. The post break slope is steeper than pre- 
dicted by a factor ~ 1/4. Moreover as 9 increases, the break 
becomes smoother (see also Tab. HI. and the time of the turnover 
increases. See text for discussion. 
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Figure 7. Same as Fig. but for a SE jet, with a comoving ve- 
locity given by Eq. |5] The jet and shock parameters are the same 
as Fig. |n] Here the post break slope given by the standard calcu- 
lation is ct2 = —p = 2.5, but again we find &2 steeper by a factor 
of ~ 1/4, immediately after the break. The lower panel shows 
that the time break increases with 9 , as well as the smoothness 
of the break, though only mildly (see also Tab. 0. See text for 
discussion. 



3.1.2 



<:>i 



Similarly the pre-break slope presents deviations from the 
standard calculations for a spherical symmetrical blastwave 
(qi = ~( 3 p~ 2 ) : Fig. |S|and Fig. [TJ lower panels, horizontal 
line on the left). If the jet is very narrow (less than few 
degrees) the lightcurve will not exhibit (depending on To) 
this first branch. For wider jets (as in our examples) the 
value ~( 3 p~ 2 ) is strictly followed only in the case of a NSE 
jet for observers around the line of sight. 



steeper. This can be understood by taking into account the 
effect of EATS. When r/T < rdj et the visible area can be 
schematized as a bright ring with radius r/T and width A 
(20% — 40% of the whole area) that surrounds a dimmer 
uniformly emitting surface (e.g. Granot et al. 1999). When 
r/T > r 9 j e t + A only the dimmer emitting surface is visible. 
As a consequence after the break the deficit in the observed 
flux is bigger and Q2 is steeper than considering, at any 
time, a uniform emitting surface. The effect is more pro- 
nounced for higher values of p since the surface brightness 
contrast between the centre and the ring increases. When 
r/T ^ rQ^t + A the lightcurve decay index should tend 
to the standard slope, because the jet surface emits almost 
homogeneously; in fact only in the case of a very narrow 
jet (6 j e t < 1°) this asymptotic slope is reached: wider jets 
lightcurves break later and they enter the trans-relativistic 
phase soon after the break, tending towards 03. 



3.1.3 a s 

On the other hand the non-relativistic regime offers a 
lightcurve branch slope independent of the spreading and 
the viewing angle. This holds for any frequency band. 
However, uncertainties in the dynamics during the trans- 
relativistic phase (see §2.1) hamper the possibility to derive 
the electron distribution power-law index even at late stages. 



3.2 Lightcurve: the time of the break 

The lower panels of Fig.|S|and Fig.[7|show, in addition, how 
the jet-break and the temporal behaviour around it change 
with the off-axis angle. As the viewing angle increases the 
transition between ai and Q2 is smoother (GL99) and the 
break time (when a ~ i2l±22i) j s retarded. 

To obtain more quantitative results we fit the 
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lightcurves F(t) with a smoothly joined broken power law 
(SBP), 



F(t) 



2 F b 



(CM*)' 



l/s 



(30) 



(Beuermann et al. 1999), where s is the smoothness param- 
eters and Fb the normalisation (Fb = F(tb) for s = 1). The 
smoothness parameters is a measure of the shape of the 
lightcurve around the break: the lower its value the smoother 
the two asymptotic slopes are joined together over the break. 
The fit is performed over the range 1 0.01 tb < t < 100 if, and 
assigning an uncertainty of 10% to each point. 

The results are summarised in Tab. Q where tb and s 
are given as a function of 8 /9j et . For NSE jet lightcurves, 
the fit gives break times that range from 1.86 days for a 
on-axis observer to 5.7 days for 9 /9j e t = 0.8, while for a 
SE jet t b = 3.57 days for 9 a = and t b = 7.44 days for 
9 /8jet = 0.8. However changing the time interval of the fit- 
ting can result in a variation of the estimated tb of a factor of 
2, while the positive correlation between tb and o generally 
holds. Tab. shows also that s decreases for larger 9 and 
the effect is greater when lateral expansion is not dynami- 
cally important (see also Fig. 4 in GL99). The bottom line of 
this discussion is that the break time is an ill-defined quan- 
tity, since the function generally used for data fits is only 
a rough approximation to the real shape of the afterglow 
lightcurve. Systematic uncertainties on the measure should 
be considered, especially when the break time is used to in- 
fer the opening angle of the jet to derive the beam-corrected 
energy output of GRBs (e.g. Berger et al. 2003). 




Figure 8. Polarisation curves for a homogeneous NSE jet, cor- 
responding to the lightcurves shown in Fig. |H] The x axis is the 
time since the trigger divided by the on-axis lightcurve time break 
tb = 1.86 days. These curves may have to be rescaled by a factor 
< 1, since Pq = 60% is taken arbitrary. For P < the vector is 
in the plane containing the line of sight and the jet axis, while 
for P > is rotated of 90° . The black solid lines correspond to 
8 = 8j et and 8 = 0.28j et . GRB afterglows always show two 
peaks, with the second one higher than the first, and with po- 
larisation angles rotated by 90°. Orphan afterglows' polarisation 
curves have two peaks, with the same angle, for 8 < 38j et that 
eventually merge into a single maximum for larger viewing angles; 
the polarisation at the peak grows with 8 a . The peak polarisation 
for an orphan afterglow can therefore be a factor ~ 2.5 larger 
than what it is expected at 8 = 8j e t- 



3.3 Polarisation curves 

The polarisation curves we obtain are given in Fig. Eland 
Fig.|i 

3.3.1 GRB afterglows 

Homogeneous GRB afterglow polarisation curves (9 /9j e t < 
1) always show two peaks, with the second higher then the 
first one, with pitch angles rotated by 90°. For 9 /9j et ~ 0.6 
this change in the direction of the polarisation vector (P — 
0) happens before the break-time measured by an on-axis 
observer, while for smaller angles it happens slightly after. 

3.3.2 Orphan afterglows 

The polarisation curve of orphan afterglows (9 /9j et > 1) 
has two peaks, with the same position angle, that eventually 
merge in a single maximum. For a NSE jet the polarisation 
at the peak grows with 9 , tending towards Po (Fig.|HJ, while 
with sideways expansion the peak value reaches a maximum 
around 9 a ~ 79j e t and then it slowly decreases; in both cases 
the polarisation peak for an orphan afterglow can be a factor 
~ 2.5 — 2.6 larger than what it is expected at 9 a — 9j e t- 

1 The break time tb is found iteratively by adjusting the fitting 
interval to the one specified and re-performing the fit until con- 
vergence. 




: 0.8fl jcl 068,., 

~io E i i i i i : 

10" 1 1 10 10 a 10 3 

t/tb 



Figure 9. Same as Fig. [S] but for a homogeneous SE jet. The 
corresponding lightcurves are shown in Fig. The x axis is the 
time since the trigger divided by the on-axis lightcurve time break 
tb = 3.57 days. All the main features are the same as for the 
NSE jet but the expected polarisation values are lower and the 
two maxima in the orphan afterglows curve merge for 8 ~ G8j e t, 
where the maximum value of polarisation is reached. Then the 
polarisation at the peak decreases. The polarisation for an orphan 
afterglow can therefore be a factor ~ 2.6 larger than what it is 
expected at 8 = 8j e t- 
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1.79 
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0.55 
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3.07 


0.60 
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Table 1. This table summarises some results from the modelling 
with a SBP (Eq.EOJ of the lightcurves shown in Fig. ©(NSE jet) 
and Fig. 171 (SE jet). The break time t\y is given as a multiple of 
the tj, fitted by an on-axis observer £(,(0)]. Note that the time 
at which the curve changes slope is postponed as the angular 
distance from the jet axis increases. The smoothness parameter 
s of the break in a lightcurve is a measure of the break shape: 
the higher its value, the sharper is the transition between the 
two asymptotic slopes. The table shows that s decreases with the 
off-axis angle. This effect is more evident for a NSE jet. 



3.3.3 Comparison with previous results 

The polarisation curves for a NSE jet with 8 < 8j e t have 
been previously published by GL99, which did not consider 
EATS. Their curves have our very same temporal behaviour, 
but their polarisation peaks are lower than ours by a factor 
from 2 to 4, depending on the viewing angle. This is the 
result of adding EATS to the computation. In this case the 
received intensity, at any given time, peaks at an angular 
distance 1/r from the line of sight, just where the linear po- 
larisation is maximised (see Eg. 1221 Eg. 1261 and Eg. 1251 1. As 
a result, the total expected polarisation is higher than for a 
homogeneously emitting surface. S99 and Granot & Konigl 
(2003) have instead explored the polarisation for a spread- 
ing jet. The first author uses a simplified model in which 
the opening angle does not change until 4 < 8j e t and then 
8 jet increases as 1/r; as a consequence he expects a third 
polarisation peak to appear at later times for large off-axis 
angles, while for small viewing angles only one peak should 
be visible. With our complete calculation of the evolution of 
the opening angle of the jet we do not obtain either of these 
effects (see Fig.|UJ. The reason is that the visible area (1/r) 
crosses both the nearest and the farthest edge of the jet 
for any off-axis angle and then 1/r remains always greater 
than 8 j et . Therefore for a HJ, afterglow polarisation curves 
always have only two peaks with orthogonal polarisation an- 
gles, for all the sideways expansion models considered in this 
paper. Granot et al. (2002) have extended the computation 
also to orphan afterglows. They obtain 

i) higher value of polarisation for GRB afterglows com- 
pared to ours with the first peak always greater then the 
second one (unlike what our curves show); 

ii) for 8 /8j et = 0.25 only one peak is visible; 

iii) an orphan afterglow peak can have a polarisation degree 
which is larger than what observed by an on-axis observer, 
but only by a factor 2. 

These effects are due to an error in their program, pointed 
out recently by the authors themselves (see Granot & Konigl 
2003). The results of their corrected code are in general 
agreement with ours (Granot private communication). 




10" 1 1 10 10 a 10 3 

t/t b 



Figure 10. Lightcurves (upper panel) and polarisation curves 
(lower panel) for a NSE jet seen at O = 0.69j e t in three differ- 
ent spectral ranges: F(v) oc i/ 1 / 3 (dot-dashed line, usually in the 
radio) F(v) oc t> - (P -1 )/ 2 (usually in the optical, solid line) and 
F(u) oc v~~ p / 2 (usually in the X-rays, dashed line). The jet and 
shock parameters are the same as Fig. |§] tt, = 3.66 days. The 
X-ray lightcurves has been rescaled by a factor 10 3 for clarity. 

3.3.4 Multiwaveband polarisation 

Polarisation due to synchrotron emission is in principle 
present at any wavelength. In Fig. 1101 we show the light 
and polarisation curves for the same GRB afterglow ob- 
served at a frequency v in three different spectral branches: 
v a <v < v m , v m <v <v a and u c < v , for an observer at 

=0.6. In the following we refer to them as the "radio" 
branch, the "optical" branch and the A"-ray branch respec- 
tively, since each waveband usually stays on that particular 
branch for most of the afterglow evolution (depending on 
fireball and shock parameters). 

Fig. 1101 shows that while the polarisation curves in the 
R and X-ray bands are very similar, the radio curve has 
a significantly lower degree of polarisation. Disticnt from 
the other lightcurves, the radio flux increases with time be- 
fore the break (F v oc t 1 ' 2 ). It should fall asymptotically as 
F v oc t -1//3 afterwards but the trans-relativistic phase sets in 
and eventually the flux rises again (F v oc t 8 ^ 5 ) following the 
non-relativistic temporal slope (Frail, Waxman & Kulkarni 
2000). 

We note that generally a spectral transition occurs in 
the radio band (8.5 GHz) at late times, when the peak 
frequency becomes smaller than the observed radio band 
[v a < v m < u ). This causes a third peak to appear in 
the radio polarisation curve, since in the radio curve joins 
the optical one; on the other hand the radio lightcurve un- 
dertakes a spectral break and eventually follows the optical 
lightcurve temporal decay (F v oc t -3 ( p-1 )' 4 before the jet 
break and F v oc t~ 3p/4 after). In our calculations for the 
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radio polarisation we do not take into account the effects 
of Faraday rotation: as long as T > 2 this effect should be 
negligible, since the electrons are all highly relativistic, with 
7 m > 100. In the trans-relativistic phase its effect can be 
sizable and it could even wash out the intrinsic small degree 
of radio polarisation. 

The X-ray polarisation curve can be affected by the 
contribution of the inverse Compton flux; in particular after 
10 — 15 days it may depolarise the synchrotron flux, depend- 
ing on the external density. However faint fluxes are usually 
observed at this epoch, too faint for detection by any X-ray 
polarimeter conceivable in the near future. In fact, detection 
of polarisation in the X-ray band will probably be possible 
only within few hours from the trigger, when the synchrotron 
flux positively dominates the inverse Compton emission. 

3.4 Summary 

In this section we have described separately lightcurves and 
polarisation curves for homogeneous jets. The variation of 
viewing angle however affects both quantities. As the off-axis 
angle increases the break time increases and since the min- 
imum in the polarisation curves occurs around tb(0 o ), it is 
delayed as well. Moreover the degree of polarisation increases 
at the peaks while the break shape becomes smoother. These 
joint characteristics should, in principle, be observable and 
be helpful for testing the model. 



4 STRUCTURED JET: RESULTS 

As a general result, our more sophisticated simulations con- 
firm all the features we described in RLR02 for the lightcurve 
of a SJ: it is very similar to the lightcurve of an HJ seen on- 
axis with same energy per solid angle and 9j e t = 2# . On 
the other hand, we show that the polarisation curves for a 
SJ present key-features that allow us to spot the underlying 
jet structure. 

4.1 Lightcurves 

The temporal index a of a SJ lightcurve is shown in the lower 
panels of Fig. [TT] and Fig. IT31 for a e =2, (3 t =l and 0v=O. 
Similar to the HJ, a SJ's lightcurve for O > few degrees 
can be approximated as a broken power-law; moreover the 
same asymptotic slopes (ay,, aa, Q3) for an homogeneous jet 
(see §3.1) describe as well the temporal behaviour of a SJ's 
lightcurve. 

It can be noticed that for viewing angles within the core 
(80 < 9 C ) the observed lightcurves are similar to those ex- 
pected from very narrow HJs: the break happens so soon 
that the first power-law branch is missing while after the 
break there is enough time to reach 02 and follow that 
slope until the non-relativistic transition sets in. For larger 
viewing angles the situation is quite the reverse: as the 9 
increases the lightcurves have time to follow more strictly 
the standard value for ct\ but there is less and less time to 
reach Q2 before entering the sub-relativistic phase. This is 
the same behaviour observed in HJs lightcurve as the jet 
opening angle increases. Another point of similarity with 
HJ's lightcurves is that (as discussed in §3.1.1 for HJs) the 
flux falls off more rapidly after the break than expected by 



standard calculations: this effects increases with the off-axis 
angle. Distinct from the HJ, the SJ lightcurves also present 
a flattening before the break that increases with the off-axis 
viewing angle. The deviations from the standard slopes be- 
fore and after the break are of the same order. 

We can conclude that HJ's and SJ's lightcurves can be 
described only roughly by broken power-laws and certainly 
only for a specific range of viewing angles. 

4-1.1 The shape of the break 

To quantify the sharpness of the break we again fit the simu- 
lated lightcurves with Eq. l30l and the results are summarised 
in Table [5] The smoothness parameter s increases and then 
saturates around 8 o /0 c = 8, where the fit yields an extremely 
sharp break (the best fit is actually with a simple broken 
power- law rather than a SBP). A NSE jet has sharper jet 
breaks compared to a SE jet, but the relation between s and 
the off-axis angle is similar. This relation is actually differ- 
ent from that discussed for the HJ (§3.2), where the break 
becomes smoother and smoother as o increases. 

4-1.2 The time of the break 

The similarity between HJ and SJ lightcurves becomes more 
quantitative if we measure the break-time as a function of 
the viewing angle. The break time t^ increases with the view- 
ing angle and the general trend can be fitted by a power law: 

t b oc (9 /9 c f 

(see Tab. 2). We underline that % oc 9 2 a is predicted by the 
structured model with a £ —2 and by the HJ model with con- 
stant total energy (once 9 D is replaced by 9j e t). 

One difference between the HJ and the SJ is that, for 
the same set of parameters, the SJ break time equals that 
of the HJ if 9jet ~ 2(9 D (RLR02) (see also Fig. ESJ. This 
factor comes from the two different origins of the break in 
the lightcurve: the HJ simply breaks when the edge of the jet 
comes into view (F ~ 9~^ t ), while the structured jet breaks 
when F c ~ O^ 1 . The following simplified calculation aims 
to show that a HJ (with the same parameters as a SJ and 
Ei SO = Ei 30 (9 )) tends to have an earlier break time if 9j et = 
9 and thus a larger opening angle (of the order of ~ 29 ) 
is needed to have a reasonable agreement between the two 
break times. The break times ratio for a non spreading jet 
is 

= f — y /3 (—) 2/ \ (31) 

tjet,h \9jet J \9 C J 

where tj etiS — Rjct - a (1 — /3cos(# D )) is the observed break 
time for SJ {Rj e t,c being the radius at which T c = l/9 ) 
and tj B t,h = - 2' jet is the break time for an HJ seen on- 
axis (Rjet being the radius at which V ~ 1/9 j e t). In Ea. 1311 
and below in Ea. 1321 we assume the same isotropic equivalent 
energy along the line-of-sight and that To does not depend 
on 9; for a spreading jet we get instead 

r^ = (/^) 8/3 ( 1+ln Sr)' (32) 

Ijet,ft \ "jet 1 \ "c I 

where the logarithmic term takes into account the exponen- 
tial behaviour of F c for F c < O^ 1 ■ If we impose then that 
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the structured and homogeneous jets lightcurves break at 
the same time Eg. I3T1 gives 1.3 < ^fi < 2.5 and Eq.glgive 



s 1.3 < ^ ^ 1.8 for 2 < |j < 40. 



4.2 Polarisation curves 



The polarisation curves of a structured jet with a E =2 and 
/3r =0 present all the same main features, regardless the 
adopted comoving lateral speed. As examples, we show a 
NSE jet f Fig. 1130 and a jet expanding with a supersonic lat- 
eral velocity Fig. 1141 . The SE jets exhibit a lower 
degree of polarisation and wider peaks than NSE jets; in 
particular for a supersonic lateral velocity, the lateral ex- 
pansion starts earlier than Ojet(r) ~ 1/r (see Fig. |5J and a 
higher degree of polarisation is present at very early times. 
Three major features characterize all the polarisation curves 
resulting from such SJs: 

• There is only one maximum in the polarisation curve. 
Since the jet is intrinsically inhomogeneous the degree of 
polarisation is greater than zero, albeit small, even at early 
times. As the visible area (~ 1/r) increases, the brighter 
inner part weights more in the computation of the total po- 
larisation (Eqs.EIOHIDand EHJ that increases until 1/r ~ 9 . 
This is the configuration with the largest degree of asymme- 
try within the visible area and therefore it coincides with a 
maximum in the polarisation curve. As 1/r becomes larger 
than 9 the degree of asymmetry decreases along with the 
polarised flux. 

• The polarisation angle does not change throughout the 
afterglow phase. Since the brightest spot is always at the 
same angle from the line of sight, the polarisation angle does 
not change through the evolution of the jet. 

• The maximum of polarisation decreases with Ei SO , The 
larger O , the larger is the visible area when 8 a — 1/r and 
the observer sees simultaneously the bright spine and the 
very dim wings so that the asymmetry is bigger. 



4-2.1 Multiwaveband polarisation 

In Fig. 1151 we compare polarisation curves (lower panel) for 
the same jet configuration (the lightcurves are shown in 
the upper panel) observed at o — 3 c in different spectral 
ranges: v a < v a < u m , v m < v a < v c and v c < v a . In the fol- 
lowing we refer to them as defined in §3.3.4. The polarisation 
in the radio branch is significantly smaller than in the opti- 
cal and X-ray bands for most of the time. As for the HJ (see 
discussion in §3.3.4), the peak frequency will eventually cross 
the radio band and the polarisation will increase and shift 
on top of the curve corresponding to v m < v < v c , while 
the lightcurve will decrease following the optical lightcurve 
slope. Different from the HJ, the optical flux has an higher 
degree of polarisation than the X-ray flux for most of the jet 
evolution. We conclude that polarisation curves depend on 
the spectrum and in particular in the radio band where two 
peaks are generally present and the degree of polarisation is 
significantly smaller than in higher frequencies bands. 
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Figure 11. Lightcurves in the R band and temporal index for a 
structured NSE jet with parameters: E c = 2 X 10 54 erg, ct e = 2, 
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eg = 0.005 and n = 1 cm -3 . The viewing angle o is indicated in 
the figure. Before the break the R band is below u c , after the break 
it is above. Thus the asymptotic slopes are: ai = — ? (p — 1) = 
-1.125, a 2 = -(3p + l)/4 = -2.125 and o 3 = -(3p - 4)/2 = 
-1.75. 
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Figure 12. The same as Fig. 1111 but for a structured SE jet 
with a comoving sideways velocity given by Eq. All the other 
parameters are the same as Fig. 1111 Here 02 = — P = —2.5. 
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Figure 13. Polarisation curves corresponding to the lightcurves 
given in Fig. 1111 The x-axis for each curve is t/(t b (^f-)), where t b 
is found by modelling the corresponding lightcurve with a SBP. 
The break times are given in Tab.|5]for 8 a > 8 C . For 8 < 8 C we 
use: t b (Q) = 0.026 days, t b (0.5) = 0.08 days and t b (l) = 1.12 days; 
these latter break times mark only the beginning of the second 
power-law branch, missing the pre-break slope. At 8 = 32 9 C the 
flattening in the lightcurve makes the measurement of the break 
time uncertain by more than a factor of 2 (see also discussion on 
the measurement of t b in §3.2) thus, despite the impression given 
by the figure, the time of the polarisation peak for 8 /8 c = 32 is 
consistent with tj,(32). 



Figure 14. The polarisation curves for a structured jet with pa- 
rameters given in Fig. 1111 that undertakes a sideway expansion 
given by Eq. llll with R = 0.1. The break times are given in Tab.l2l 
for 8 > 8 C . For 8 a < C we use: t b (0) = 0.06 days, 4(0.5) = 0.092 
days and t b (l) = 0.3 days; these latter break times mark only the 
beginning of the second power-law branch, missing the pre-break 
slope. In this case the corresponding lightcurves show a more 
pronounced flattening and the determination of the break time 
becomes more precise. All the times of the polarisation peaks are 
consistent with the break times in the lightcurves, since they are 
shifted in the figure by less than a factor of 2 (see discussion on 
the measurement of t b in §3.2). 





NSE 




SE Eq.EI 


SE Eg. [HI 


d / 8 C 


H 

*b(2) 


s 


ib(2) 


s 


h 

ib(2) 


s 


2 


1 


1.29 


1 


1.01 


1 


2.63 


4 


5.53 


3.72 


1.72 


0.68 


4.50 


4.04 


8 


31.2 


12.19 


7.66 


1.71 


22.30 


12.91 


16 


187.53 


20 


37.22 


20 


119.50 


20 


32 


846.33 


20 


171.74 


20 


699.06 


20 


64 


4602.53 


20 


739.32 


20 


3762.45 


20 



Table 2. This table summarises some results from the modelling 
with a SBP (Eq. |3UJ of the lightcurves shown in Fig. [HJ (NSE 
jet) and Fig. 1121 (SE jet) for a structured jet. The values are 
given for O > 9 C since for smaller viewing angles the lightcurve 
can only be fitted very roughly by a SBP, missing the pre-break 
branch. The break time t b increases with the viewing angle and 
the general trend of all the measured break times can be fitted by 
a power-law t b oc (8 /8 c ) 2 . The smoothness parameter s increases 
and then saturates around 8 /8 c = 8. A NSE jet has sharper jet 
breaks compared to a SE jet. 



5 GAUSSIAN JET 

We end our investigation of GRB jet luminosity struc- 
tures with a brief discussion of the Gaussian jet curves, 
which present features of both the HJ and the SJ (see also 
Salmonson et al. 2003). 



5.1 Presentation and Dynamics 

Perhaps a more realistic version of the standard "top hat" 
model is a jet whose emission does not drop sharply to zero 
outside the characteristic angular size (9 C ). Such a configu- 
ration can be described with a Gaussian distribution of the 





Figure 15. Lightcurves (upper panel) and polarisation curves 
(lower panel) for a structured NSE jet seen at 8 = 38 c in three 
different frequency ranges: v a < v < v m (dot-dashed line) u m < 
Vo < Vc (solid line) and v c < v (dashed line). All the other 
parameters are given in Fig. 1111 t b ~ 1 day. 
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Figure 16. Lightcurves in the R band and temporal index for 
a Gaussian NSE jet with parameters eo = 10 53 erg, Fq = 10 4 , 
C = 10°, e jet = 90°, e e = 0.01, e b = 0.005 and n = 10 cm" 3 . The 
viewing angle 9 is indicated in the figure in units of 8 C . In the 
upper panel, lightcurves with O < 9 C are indistinguishable and 
are therefore labelled all together. The solid curve has 9 = C . 
Horizontal dashed lines in the lower panel indicate the asymptotic 
slopes in the various regimes. 
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Figure 17. Polarisation curves corresponding to the lightcurves 
shown in Fig. 1161 The x-axis is plotted in units of t b (0), the break 
time of the on-axis lightcurve, as in Figs. and |H1 

energy per unit solid angle (Eq.EJ. In this jet the luminosity- 
varies slowly within 9 C (t^tj — 1-65) and it decreases expo- 
nentially for 9 > 9 C . In principle, the dynamics of a Gaussian 
jet could be non-trivial due to the fact that, unlike the uni- 
versal SJ (a e = 2) and, of course, the HJ matter at different 
angles does not start spreading at the same radius (in partic- 
ularly for 9 > 9 C ); this can develop lateral velocity gradients 
and transversal shock waves. For this reason we restrict our- 
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Table 3. The same as Tab. 1 for the lightcurves shown in Fig. 1161 
(NSE jet). The break time is given as a multiple of the if, fitted 
by an on-axis observer [tj, (0) = 4.56 days]. 



selves to the case of a non-spreading jet. This assumption is 
supported by recent numerical hydrodynamical simulation 
for a Gaussian jet evolution; they suggest that the trans- 
verse velocity remains below the sound speed as long as the 
evolution is relativistic (Granot & Kumar 2003). 



5.2 Lightcurves 

The resulting lightcurves are shown in Fig. 1161 It shows that 
for 9 < 9 C the lightcurves are very similar to a HJ's ones: 
the slopes, the values and the behaviour of the break time 
as a function of angle (see Tab.[3J are consistent with what 
expected for a NSE HJ seen within the cone (see Tab.0. On 
the contrary the break does not become smoother (3rd col- 
umn of Tab.[j|J as the line of sight approaches the edge, like 
in the HJ but its shape remains rather unaltered. Finally, 
on average the Gaussian jet has a less sharp break in the 
lightcurve. This is in agreement with previous calculations 
(Granot & Kumar 2003). For 9 > 29 c instead the pre- break 
slope becomes flatter and flatter and for 9 > 89 c it becomes 
positive; in this latter case the lightcurve is actually domi- 
nated by the emission coming from the core. 



5.3 Polarisation curves 

The polarisation curves present intermediate characteristics 
between the HJ and the SJ ones (Fig ITTJl. The absence of 
edges and the presence of a symmetric luminosity gradient 
with respect to the jet axis, produce (as in the case of a SJ) 
a one-peak curve with a constant polarisation angle. On the 
other hand, the exponential decrease of luminosity outside 
the core makes the relation between polarisation curve and 
lightcurve more close to the HJ's one: the peak is located 
after the break time in the total flux. As the viewing an- 
gle increases, especially for 9 > 9 C the maximum in the 
lightcurve moves torwards t = tb, but eventually the core 
starts to dominate the lightcurve since early times because 
of the exponential luminosity distribution and the polarisa- 
tion curves resembles that of an orphan afterglow. 
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6 COMPARISON AND DISCUSSION 

The previous sections confirm what first claimed in RLR02: 
the lightcurves of a HJ and of SJ with a e = 2 are very simi- 
lar but their intrinsic features depend on the viewing angle 
for the SJ and on the opening angle for the HJ. On the other 
hand the polarisation curves are completely different. These 
characteristics are the direct consequences of an energy dis- 
tribution e oc 6~ 2 : the lightcurve is dominated by the line 
of sight emission while (for any a e ) the polarisation curve is 
dominated by the emission coming from an angle 1/r with 
the line of sight. This means that the total flux we receive 
does not bear footprints of the jet structure while the ob- 
served polarisation does. In Fig. 1181 we directly compare the 
lightcurves and the polarisation curves from HJ and SJ with 
Ei so (8 ) = Ei 30 and 6j e ± = 20 o . These are the parameters 
for which their lightcurves are more similar. We also show 
the lightcurve and polarisation curve for a Gaussian jet with 
Ei ao (9 ) — Ei SO and 8 C ~ 0.6#j e t. 



6.1 Lightcurves 

Fig. ED (upper panel) summarises the comparison discussed 
§4.1 and §5.2 among the characteristics of the lightcurves of 
the three jet structures. In particular it should be noticed the 
flattening before the break, present in the SJ lightcurve and 
the almost perfect match between the HJ and GJ lightcurves 
(the GJ lightcurve has been divided by a factor 2 or it would 
be overlaid don the HJ one). The pre-break bump is actually 
the only sign of the underlying jet structure: when 1/r ~ O 
and the lightcurve breaks, the jet core is also visible and 
its contribution gives that flux excess that we perceive as 
a flattening. The larger is the viewing angle the smaller is 
Ei SO and more prominently the core out-shines the line of 
sight flux at the break. Some authors claimed (e.g. Granot 
& Kumar 2003) that this feature can used to discriminate 
between the HJ and the SJ. However the shape and the in- 
tensity of the flattening depends on how the wings and the 
core join together and this is a free parameter (/3 e ). Con- 
sequently comparing the model with observations provides 
a way to fix the shape of the energy distribution but not a 
way to test the model itself. Besides the pre-break flattening 
the temporal behaviour of the lightcurves plotted in Fig. 1181 
is identical and for this reason the same data can be fitted 
with any of the three models (e.g Panaitescu & Kumar 2003, 
for SJs and HJs). 



6.2 Polarisation curves 

The comparison in Fig. 1181 (lower panel) clearly shows the 
main differences between the polarisation curves of a SJ, of 
a GJ and of a HJ: 



• Different behaviour from the beginning. 

The HJ and GJ do not display polarisation at early times 
but only when the visible area intersects the edge of the jet 
and an asymmetry is present. However, the SJ shows, from 
the beginning, regions with different luminosity within 1/r 
from the line of sight and therefore we observe a non zero, 
albeit small, degree of polarisation even at very early times. 

• Different evolution of the polarisation angle: 90° rota- 
tion vs constant 
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Figure 18. Lightcurves (upper panel) and polarisation curves 
(lower panel) comparison between a SJ seen at O = 4° (the other 
parameters are given in Fig. 1111 . a homogeneous jet with 8j e t ~ 
29 , Ei so ~ Ei ao (9 ) and To ~ rrj(t9 ) and a Gaussian jet with 
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The polarisation angle for an HJ rotates by 90° at ~ ti\ 
this is because opposite parts of the jet dominate the total 
polarisation when 1/r overtakes the nearer edge and when 
it reaches the furtherest one. For a jet seen off-axis the break 
in the lightcurve is not sharp (§3.1.1) and the time break if, 
occurs roughly when r ~ (0 o + 9jet)/2, midway between the 
two edges. That is why the change in the polarisation angle 
corresponds to the jet break. This behaviour is observed for 
any off-axis angle, for all sideways expansions considered 
in this paper. On the other hand in the SJ and GJ the 
same spot (towards the core) always dominates the total 
polarisation and consequently the angle remains constant. 

• Different number of peaks: 2 vs 1 
The HJ has two peaks in the polarisation lightcurve, with 
the second always higher than the first. This is of course 
related to the presence of the two edges of the jet that are 
the source of the asymmetry within the visible area (Fig. 
1 and 2 GL99). The SJ has only one maximum when the 
observer sees light coming from the core at 1/r from the 
line of sight. This is when a break occurs in the lightcurve 
and the polarised flux is almost completely dominated by the 
flux coming from the more powerful region of the jet. The 
behaviour of the GJ is in between: the polarisation angle is 
constant as for a SJ, but the polarisation peak is delayed 
with respect to the break time, as in a HJ. 
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7 CONCLUSIONS 

We have presented a thorough analysis of light and polari- 
sation curves from GRB jets. Our study is based on codes 
that numerically integrate the jet equations taking into ac- 
count the equal arrival time surfaces and the finite width of 
the radiating shell. We initially concentrate on homogeneous 
jets, pointing out some features of the jet evolution which 
had been overlooked in previous works. In particular we un- 
derline the difficulty of measuring the break time and the 
electron energy distribution index p modelling the lightcurve 
with a SBP; we also point out that in the non-relativistic 
regime the usually assumed Sedov- Taylor solution could be 
not correct when applied to narrow jets. Concerning polar- 
isation, we find the presence of only two peaks in the po- 
larisation curve, for any off-axis angle and for any sideway 
expansion velocity considered so far in literature. 

The main result of this paper is that the polarisation 
temporal behaviour is found to be very sensitive to the lu- 
minosity distribution of the outflow unlike the total flux 
curve. We have achieved the goal of calculating, for the 
first time, the polarisation from a universal structure jet and 
from a Gaussian jet and performing a comparison (of both 
lightcurves and polarisation curves) with what is predicted 
by the standard jet model. 

In particular the derivation of the light and polarisation 
curves for structured jets has been performed under several 
assumption for the dynamics. For simplicity we study only 
a non-spreading GJ. Our numerical lightcurves confirm the 
previous conclusion (RLR02; Granot & Kumar 2003; Ku- 
mar & Granot 2003; Salmonson 2003) that, based on the 
lightcurve properties, it is extremely hard to infer the struc- 
ture of the jet. Polarisation curves, on the other hand, are 
extremely different. Since the brightest part of the jet is al- 
ways on the same side for structured (SJ and GJ) outflows, 
the position angle of polarisation remains constant through- 
out the whole evolution. In addition, for a SJ the polarisa- 
tion peaks coincident in time with the jet break in the light 
curve, which instead corresponds to the time of minimum 
polarisation in homogeneous jets. The exponential wings in 
a GJ, however, shift the position of the peak after the break 
in the lightcurve, a feature that marks the difference be- 
tween the GJ and the SJ predicted polarisation. We should 
however stress that due to the many uncertainties inherent 
in the derivation of polarisation curves, it is hardly possible 
to use them to measure in a fine way the energy distribu- 
tion of the jet (e.g. tell a e oc 6~ 2 structured jet from a 
e oc #~ 2 ' 5 one). What polarisation can robustly determine 
is whether the energy distribution in the jet is uniform or 
centrally concentrated. Alternative approaches, such as the 
observed luminosity function (RLR02) are also important to 
further constrain the jet profile, even though the data seem 
not to be accurate enough at this stage (Perna et al. 2003; 
Nakar et al. 2004). 

Finally we present unprecedented polarisation curves 
corresponding to different spectral branches, both for a HJ 
and for a SJ: we find a spectral dependence of the degree 
of polarisation and we expect changes in the temporal be- 
haviour of the polarised fluxes (in all bands but in particu- 
lar in the radio band) associated with spectral breaks in the 
lightcurve. 

It should be emphasised that these results hold in the 



absence of inhomogeneities in the external medium and/or 
in the luminosity distribution within the jet (other than 
e oc 9 a ' ). Any breaking in the fireball symmetry causes ran- 
dom fluctuations of both the value and the angle of the 
polarisation vector (Granot & Konigl 2003; Lazzati et al. 
2003; Nakar & Oren 2004). In the lightcurve the main con- 
sequence is the presence of bumps and wiggles on top of the 
regular power-law decay (e.g. Lazzati et al. 2002; Nakar, Pi- 
ran & Granot 2003). Actually a complex behaviour for the 
lightcurve has always been observed so far together with a 
peculiar polarisation curve, such as in GRB 021004 (Rol et 
al. 2003, Lazzati et al. 2003) and in GRB 030329 (Greiner 
et al. 2003). On the other hand GRB 020813 presents an ex- 
tremely smooth lightcurve (Gorosabel et al. 2004; Laursen 
& Stanek 2003) associated with a very well sampled polari- 
sation curve (Gorosabel et al. 2004) that is characterised by 
a constant position angle and a smoothly decreasing degree 
of polarisation. This burst is thus particularly suited for a 
proper comparison with data of the models described in this 
paper. Lazzati et al. (2004) has performed the modelling of 
the polarisation curve according to several models (includ- 
ing all the models considered here) and they find that the 
structured model can successfully reproduce the data and 
it can predict the jet-break time in agreement with what 
measured in the lightcurve. 

Further complication can arise from the presence of a 
second non negligible coherent component of the magnetic 
field in the ISM. Our results have been obtained assuming 
that the magnetic field responsible for the observed syn- 
chrotron emission is the one generated at the shock (thus 
tangled at small scales). This is a reasonable assumption 
since the compression of a standard interstellar field is far 
too low to produce the observed radiation. However if the 
burst explodes in a pre-magnetised environment or if the 
jet is magnetic dominated and the field advected from the 
source survives till the afterglow phase, then the polarisation 
curve will be the result of the relative strength of the two 
components of the magnetic field. Granot & Konigl (2003) 
have discussed the polarisation curve in the former case for 
an homogeneous jet propagating through a magnetic wind 
bubble. 

Finally the intrinsic polarisation curve of the afterglow 
can be affected by the dust present both in the Milky Way 
and in the host galaxy. Lazzati et al. 2003 discuss the mod- 
ification of the transmitted polarised vector and they show 
that in GRB 021004 a sizable fraction of the observed po- 
larised flux is likely due to Galactic selective extinction. 

Despite the difficulties inherent to polarisation observa- 
tions and modelling (Lazzati et al. 2004), we believe that 
polarimetric studies are of great importance in determining 
the structure of GRB outflows. 
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